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Abstract. 

The analytical and numerical solutions of the Mathisson-Papapetrou equations 
under the Mathisson-Pirani supplementary condition describing highly relativistic 
(ultrarelativistic) motions of a spinning particle in a Schwarzschild field arc 
investigated. The known condition \So\/mr <C 1, which is necessary for a test particle, 
holds on all these solutions. The explicit expressions for the non-equatorial circular 
orbits, in particular for the space boundaries of the region of existence of these orbits, 
are obtained. The dynamics of the deviation of a spinning particle from the equatorial 
ultrarelativistic circular orbit with r = 3M caused by the non-zero initial value of the 
radial particle's velocity is studied. It is shown in the concrete cases that spin can 
considerable influence the shape of an ultrarelativistic trajectory, as compare to the 
corresponding geodesic trajectory, for the short time, less than the time of one or two 
revolutions of a spinning particle around a Schwarzschild mass. 



PACS numbers: 04.20.-q, 95.30.Sf 
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1. Introduction 

The present paper deals with the Mathisson-Papapetrou equations describing the motion 
of a spinning test particle in a gravitational fieldf |2]. Partial solutions of these 
equations have been considered since 1951 |4J. Different aspects of the orbital motions of 
a spinning test particle in the gravitational fields were investigated in many papers [5-41]. 
For the correct description of the trajectories of a spinning test particle by the Mathisson- 
Papapetrou equations it is necessary to take into account the condition \So\/Mm <C 1 
(| So I and m are the absolute value of spin and the mass of a particle respectively, and M 
is the mass of a gravitational source). Probably, R.Wald was the first who introduced 
this condition explicitly Another point of importance concerning the Mathisson- 
Papapetrou equations is the supplementary condition for the concretization of the center 
mass position of a spinning particle. In relativity, the position of the center of mass 
of a rotating body depends on the frame of reference [42-44]. CMoller distinguished 
the proper center of mass and the non-proper centers of mass. Because the correct 
definition of the center of mass for a spinning particle is a subject of discussion, different 
conditions are used. Most known are the Mathisson-Pirani pfl ES] and the Tulczyjew- 
Dixon conditions [JOl Elj (the first of them is often called as the Pirani condition, though 
before [15] it was used in [I]). More rare is the Corinaldesi-Papapetrou condition 
(All these conditions are described in section 2). For example, the Mathisson-Pirani 
condition was used in [8, 19, 33, 40] and the Tulczyjew-Dixon condition was used in [9, 
12, 24-26, 28-30, 32, 35, 36, 39, 41]. Authors of some papers take into account both 
the Mathisson-Pirani and Tulczyjew-Dixon conditions [18, 20, 37, 38]. More about the 
supplementary conditions we shall write below. Here we stress that it is important to 
study the solutions of the Mathisson-Papapetrou equations under different conditions 
and to compare their physical consequences. 

The feature of importance of the Mathisson-Papapetrou equations is that they are 
obtained without the restriction on the particle's velocity. This point is common for 
these equations and the geodesic equations. However, it is easy to see the essential 
difference: the highly relativistic orbits of a spinless test particle following from the 
geodesic equations have been studied in full detail in different gravitational fields whereas 
the possible orbits of a fast particle following from the Mathisson-Papapetrou equations 
are investigated insufficiently even in a Schwarzschild field. The reason for this situation 
is connected with the a priori assumption that the world lines of a spinning test particle 
are close to the corresponding geodesic lines in the whole range of the particle's velocity. 
Such an assumption is acceptable for the low velocity of a spinning particle but is not 
evident for the high velocity when a priori we cannot exclude the specific influence of 

f It was pointed out in some papers that these equations were first derived by A.Papapetrou. Perhaps, 
such a conclusion was inspired by the non-adequate sentence from [2], p. 254. About the role of Dr. 
Myron Mathisson (14.12.1897, Warsaw - 13.09.1940, England; P.Dirac published Obituary in Nature 
1940, 146, 613) and his dramatic life see After pQ these equations were rederived and reformulated 
in more than 20 papers, many of them are cited in ^5] [28] 
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the great relativistic Lorentz 7- factor on the interaction of the spin with a gravitational 
field. 

It has been shown in [SUES! that, from the point of view of the observer comoving 
with a fast spinning test particle, the 3-acceleration of this particle relative to the 
corresponding spinless particle (i.e. the particle with the same initial values of the 
coordinates and velocity) in a Schwarzschild field is proportional to r y 2 . In this sense 
we say that the highly relativistic motions (with r y 2 ^> 1) of a spinning test particle 
essentially differ from the geodesic motions. (Other possible features of the "essentially 
non-geodesic motions" of a spinning particle we shall consider below). The force 
deviating the motion of a spinning particle from the geodesic motion is the gravitational 
ultrarelativistic spin-orbit force [23 E3] • 

Naturally, further investigations must be carried out to determine the possible 
deviation of the orbital motion of a spinning particle from the corresponding geodesic 
orbits (i.e. the orbits with the same initial values of the coordinates and velocity) from 
the point of view of the observer at rest relative to a Schwarzschild source. 

The purpose of this paper is to present the new partial solutions of the Mathisson- 
Papapetrou equations under the Mathisson-Pirani supplementary condition in a 
Schwarzschild field describing the highly relativistic motions of a spinning test particle. 
We shall consider both the analytical and numerical solutions. In this context we 
stress that the wide classes of the numerical solutions of the Mathisson-Papapetrou 
equations under the Tulczyjew-Dixon condition in the Schwarzschild and Kerr fields were 
investigated in [24-26, 30, 35, 39]. The results of these papers are of importance for study 
of the dynamics of the binary black holes. We point out some differences between these 
results and our approach: (1) we use the Mathisson-Pirani supplementary condition and 
argue this choice; (2) many results of [24-26, 30, 35, 39] describe the situations when the 
value \So\/Mm is of order 1, i.e. when the Wald condition is not satisfied. Whereas all 
our solutions satisfy the relationship \So\/Mm < 1; (3) we focus our attention on the 
highly relativistic motion just of a microscopic spinning particle, because in reality any 
macroscopic test body cannot be considered as highly relativistic. A common feature 
of all our solutions is that they describe the ultrarelativistic motions with the 7- factor 
satisfying the relationship: 7 -2 is of order | Sq\ /Mm. 

This paper is organized as follows. In section 2 we review the Mathisson-Papapetrou 
equations and the supplementary conditions. The partial analytical solutions of these 
equations describing non-equatorial circular orbits of a spinning test particle in a 
Schwarzschild field are studied in section 3. In section 4, using the integrals of the 
energy and angular momentum in a Schwarzschild field, we write the strict Mathisson- 
Papapetrou equations under the Mathisson-Pirani supplementary condition in the form 
of the second-order differential equations relative to the coordinates r and tp for the 
equatorial motions of a spinning particle. The numerical solutions of these equations 
are investigated in section 5. 

Throughout this paper we use units c = G = 1. Greek indices run 1,2,3,4 and latin 
indices 1,2,3; the signature of the metric (-,-,-,+) is chosen. 
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2. The Mathisson-Papapetrou equations and supplementary conditions 

The Mathisson-Papapetrou equations can be written in the form [IJ Ej 

up ds~ =0 

ds ds ds 

where u x is the 4- velocity of a spinning particle, S^ u is the antisymmetric tensor of spin, 
m and D/ds are, respectively, the mass and the covariant derivative with respect to 
proper time s; R^ pcr is the Riemann curvature tensor of the spacetime. The Mathisson- 
Pirani supplementary condition for equations (JTJ, (J2J) is [T| 143] 

S^Uu = 0. (3) 

(Relationship (j3J) was used earlier in special- relativistic electrodynamics [3%]). 
The Tulczyjew-Dixon supplementary condition is |"""""J Ej 

S^P„ = (4) 

where 

DS UfJl 

P» = mu » + U (5) 
ds 

is the particle's 4-momentum. 

Papapetrou and Corinaldesi used the condition jlj 

S i4 = 0. (6) 

(It is easy to see that for m = relationship (0) coincides with (jSJ)). 

It is known from [JOl EHl that in the Minkowski spacetime the Mathisson- 
Papapetrou equations under condition have, in addition to the usual solutions 
describing the straight worldlines, a family of solutions describing the helical worldlines. 
(As a partial case, this family contain the circular solutions). These unusual solutions 
often are called as the Weyssenhoff orbits. By CMoller, the usual solutions describe 
the motion of the proper center of mass of a spinning body (particle), and the helical 
solutions describe the motions of the family of the non-proper centers of mass [4*3*1 144j . 

To avoid the superfluous solutions of equations ((H), (j2J), instead of condition (jHJ) the 
Tulczyjew-Dixon condition was introduced. This condition picks out a unique worldline. 
However, the question arises: is this worldline close, in the certain sense, to the usual 
worldline of equations (J""}, (J2J) under condition (jHJ)? It is simple to answer this question 
when the relationship 

DS V » 



m\u v \ 3> 



(7) 



takes place, because in this case conditions (JU) practically coincides with (jSJ). It is easy 
to check that relationship ((7|) holds if the Wald condition 

|e|<l, e = -^- (8) 
1 1 Mm K ' 
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is taken into account and, in addition, if the particle's velocity is not too close to the 
velocity of lightf . Here we point out that there is the connection between Sq snd S^ v 

Si = \s»„S^. (9) 

Another situation with relationship (JJJ) takes place for the highly relativistic motions 
of a spinning particle when u v is proportional to the 7-factor, because the term DS ull /ds 
in (J7|) depends on u v and on 7-factor as well. Therefore, in general, for the ultrtelativistic 
particle's velocity it is not obvious that (J7j) is valid. 

So, if the velocity of a spinning particle is not too high, one can "forget" about the 
Mathisson-Pirani supplementary condition, because it is sufficient to use the Tulczyjew- 
Dixon condition. Whereas for highly relativistic motions the Mathisson-Pirani condition 
must be taken into account. We stress that just the Mathisson-Pirani condition is 
derived in some papers by different method [HI E2 EH] • (For example, we agree with 
the statement that the Mathisson-Pirani condition "... arises in a natural fashion in 
the course of the derivation", [22], p. 112). That is, the Mathisson-Pirani condition 
is the necessary one, though often, with high accuracy, it can be substituted by the 
Tulczyjew-Dixon condition. 

We summarize: the existence of the superfluous solutions of the Mathisson- 
Papapetrou equations under the Mathisson-Pirani supplementary condition is not a 
reason to ignore this condition. Because just among all solutions of these equations 
under Mathisson-Pirani condition there is the single solution describing the motion of 
the particle's proper center of mass. As we shall see in section 5, it is not a great problem 
to recognize such a solution, at least in some cases of importance. 

In this paper we shall use the Mathisson-Pirani condition. At the same time, in 
some points, we take into account condition (J3J). More exactly, for the concrete solutions 
of the Mathisson-Papapetrou equations we shall estimate the value of the left-hand side 
of relationship (J3J) . 

3. Non-equatorial circular orbits of a spinning particle in a Schwarzschild 
field 

The subject of this section is inspired by the discussion in papers [2H E2] where the 
question was considered concerning the existence or non-existence of non-equatorial 
circular geodesies with constant latitude in a Kerr metric. It was shown in [22] that 
such geodesies do not exist. An interesting question is the following: are there not the 
non-equatorial circular orbits with constant latitude in the Schwarzschild or Kerr metric 

f For different estimates it is heipful to write numbers of |e| for various particles. It is known that 
in units where c = G = 1 the numerical values of the electron's mass, the Sun's mass and the Planck 
constant are equal respectively to: m = 2.3 • 10 -66 , M = 5- 10~ 6 , h = 2 ■ 10~ 86 (see, e.g. [HI]). Then for 
a Schwarzschild black hole of mass that is equal to three of the Sun's mass we have |e e | = 4.6 • 10 -17 . 
The analogous values for a proton and a neutrino (of mass corresponding to the energy of order leV) 
are: |e p | = 2.5 • 10~ 20 , \e v \ — 2.3 • 10 _n . By these values it is easy to recalculate the corresponding 
values for a massive Schwarzschild black hole or for a hypothetic microscopic black hole. 
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according to the Mathisson-Papapetrou equations? Now we shall consider the case of a 
Schwarzschild metric. 

Here and in the following the standard Schwarzschild coordinates x 1 = r, x 2 = 
9, x 3 = ip, x A = t are used, and we put r > 2M (M is the Schwarzschild mass), i.e. 
the region above the horizon surface is under investigation. 

Let us check, by the direct calculations, do equations (JTJ), (J2J) have the solutions 

with 

r = const ^ 0, 6 = const ^ 0, ir/2, ir, (10) 

o dip a dt . , . 

u 6 = -f- = const ^ 0, u = —— = const ^ 0. (11) 
ds ds 

Relationships (|T(H) . (jll|) correspond to the circular motions with the constant particle's 
velocity. 

We start from equations (0). 

3.1. The consideration of equations (0) 

It is convenient to use the representation of equations ((21) through the 3-vector of spin 
Si where by definition 
1 



Si = ^V=ge ikl S kl (12) 
{eiki is the spatial Levi-Civita symbol). The inverse relationship is 

S kl = -^e klm S m . (13) 



-g 

By condition (jHJ) we have 

(We point out that in many papers the 4-vector of spin s\ is used where 

sx = ^V=ge Xllva u' t S vv . (15) 



The relationship between Si and s\ is Si = MjS 4 — u^Si 

According to jlU] three independent equations from (J2J) may be written in the form 

u A Si - ii 4 Si + 2 (ii[iUi] - u n u p r^ [4 Ui]j S k u k + 2S n T^u i] u 7r = 0. (16) 

Here and in the following a dot denotes usual differentiation with respect to proper time 
s, and square brackets denote antisymmetrization of indices. We stress that equations 
under condition (J3J) is equivalent to equations (|TB|) . 

Taking into account equations (fTUj) . (fTTj) we write three equations from (fTB^) as 

& + 5 3 n 3 n 4 M4 (rf 4 - rf 3 ) = 0, (17) 

5 2 - s 3 wVu 4 rl3 = o, (18) 

5 3 + S lU3 (g 44 Tl 4 - g 33 Tl 3 ) - S 2 T 2 33 u 3 g 33 = 0. (19) 
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It is easy to see that the set of equations (|17j) -(|19| ) has the partial solution with 

S 3 = S v = 0, Si = S r = const, S 2 = Sq = const, (20) 

(21) 



Using the explicit expressions for g^ u and V* in the standard Schwarzschild coordinates 
we rewrite equation (|2"Tj) 



r J r sin 6 

In the following we shall analyse equations (0) for the case just when relationships (|2T)|) . 
take place. 



3.2. The consideration of equations (QJ) 

Let us write equations (JTJ) taking into account relationships (fTU|). (fTTj) . (|T3|). ((T^J) • Using 
the explicit expressions for R^ pa , two equations of set with A = 1 and A = 2 can be 
written as 



m(r^V + r^uV) + u%Tl 3 - g 44 g 33 T 



3/-pl 



,44 „ -nl 



1 



44/ 1 



:(gMT u ll u S 2 + 



g 33 Tl z u- i u- i S 2 - g^T^u-'u-'Sx 
1 



3M 



u 3 S<? sin 6* 



(23) 



-3 

3M 



u 3 S± sin 0. 



(24) 



Two other equations of set (0), with A = 3 and A = 4, become the identities. By the 
explicit expressions for T^ a we write the linear combinations of equations (J2HJ) . ((211): 



MS 



(« 3 ) 2 (— -SiCot#)sin#-mw 3 - 

r r 4 sin 9 



2 V) 2 



3u 

—m(u 3 ) 2 sin 2 # H — 

4\2 



m / 2M 
+ — I 1 I ( " 



5*2 — Sir I 1 ) tan H 



0. 



3Mg x 
r 3 cos # ' 

sin #+ 



(25) 



v / ■ (26) 

Using relationship (|2*2|l . the set of second-order algebraic equations (relative to w 3 , u 4 ) 
f)25jl . can be written as 

ai(u 3 ) 2 + a 2 u 3 + a 3 (u 4 ) 2 = 0, (27) 

b 1 (u 3 ) 2 + b 2 u 3 + b 3 (u 4 ) 2 = d (28) 



where 



6S 2 



1 - 



2i\i 
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-i 

, b 2 = -m, b 3 ■■ 



So 



3M 



rsin# \ / 
3MS 2 



sin 2 9 



3M 



r 



MSo 



r 4 sin 9 ' 



d 



r 4 sin 9 



1 - 



3M 



(29) 



In addition to equations (J27|) . (|28|) . the components of the particle's 4- velocity u 3 , w 4 
must satisfy the general relationship = 1 which in our case takes the form 

g 33 (u 3 ) 2 + g A4 (u 4 ) 2 = l. (30) 

We point out that in (|2H|) r ^ 3M because according to equations (j22J)-(I21I) at 
r = 3M we have the two possibilities: (i) cos 9 = 0, i.e. the case of the equatorial 
circular motion which was considered in [10], and (ii) 5*2 = 0, Si ^ 0; however, it is easy 
to check that in the second case the expressions for u 3 ,w 4 following from ()25|) . (|2*Bj) do 
not satisfy equation (j30|) . 

Do equations (J27j) . (J25J) . (J3T?|) have solutions at r ^ 3M? To answer this question, 
let us substitute the expression for (w 4 ) 2 through (w 3 ) 2 from (j3Uj) in (|2*7j) . (|2~8*)). Then we 
have 



/ 3\2/ #33 x . 3 i a 3 n 

(it ) (ai - a 3 — ) + a 2 u H = 0, 

944 944 

(n 3 ) 2 ^-^— ) + & 2 « 3 + — -^ = 0. 

944 944 

Taking into account ()29|) we have a\ — a 3 g 33 / g AA = 0, and it follows from (f3*Tj) 



mr 



u 



6S2 sin 9 
By P3), (EH) we obtain 



OS, ^ o). 



mr 



1 



2M 



-1/2 



1 + 



36S| 
m 2 r 4 



2\ 1/2 



(31) 
(32) 

(33) 
(34) 



Substituting ([3llj) in (EU) we get the condition under which three equations (|27|). 
(I3T?|) are compatible 



2 2 
m r 



36S 2 sin 2 



S 2 



r sin 9 



3M 



sin 2 fl 



+ 



m 2 r 



MS 2 



1 - 



1 - 



2M 



3M 
r 

-1 



-1 



MS 2 sin 9 



1 - 



2M 



-1 



6S2 sin # r 4 sin 9 \ r 
Equation (J33J) can be solved relative to sin 2 9: 



3MS 2 
r 4 sin 9 



1 - 



3M 



-1 



(35) 



sin 2 fl 



1 



2M 



' M , 
— 4 



9M 



1 + 36 



m 2 r 4 



-6 1 



2M 



1 - 



3M 



(36) 



So, equation (1) in a Schwarzschild field at r 7^ 3M has the partial solutions with 

O, (El, (ESD- 
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3.3. The spatial region of existence of these solutions 

Now we shall write the relationship between S2 and Sq. By (j9j)- fTTf . (J20|) . (}22|) we find 



C2 

~~ ^2(„A\2 



r 2 (w 4 ) 

Taking into account expression 



1+1- 



2M 



1 



3M 



cot 2 fl 



(37) 



we obtain from (|37|) 



\S n 



mr 



6SI 
m 2 r 4 



x 



-6 ^1 

2M 



2M 



3M 



36 : 



M 



3M 



1 + 36 



r 

sj 

m 2 r 4 



9M 



Si 



m 2 r A 



1/2 



-1/2 



(3? 



We remind that the condition 
(JHJ) and (|38p. it is necessary 



m 2 r 4 



<C 1. 



must be fulfilled. Therefore, according to equations 

(39) 



Then by (JHD, (123) we write 
\Sq\ 



mr 



m 2 r A 



3M 



-1/2 



(40) 



Naturally the values of the right-hand side of equation must be between and 
1. According to (J36)) this condition is fulfilled if 



15 

y 



M 1 



6 Si 



5 m 2 r 4 



< r < 3M 1 



6^ 2 2 
m 2 r 4 



By equation ()39|1 relationship (|4~T|) can be written approximately as 
15 



7 



-M <r < 3M. 



(41) 



(42) 



This relationship determines the space region where equations have the partial 

solutions for which expressions (fTU j) . (fTTjl. pljl . (j^ . pi ) . pHjl . (Pjl hold. 

According to equation (|3^j) and (|3^j) the expression for sin 2 9 can be written 
approximately as 



-6 1 



sin 2 # 
2M 



2M 



M 
r 



4 - 



9M 



3M 



(43) 



By f!43j) it is easy to find the interval of changes of sin 2 9 when r is changed in interval 
) . The minimum value of sin 2 9 is 



sin 2 fl| 



0.465 



and this value is achieved at 
5V5 



2.35M. 



(44) 



(45) 



3V5- 1 

At r = 2.35M and r = 2.5M we have from equation (|43j) the value sin 2 9 = 0.5, 9 = 45°. 
At r = 15M/7 and r = 3M the value of sin 2 9 is maximum, namely sin 2 9=1. 
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3.4- The interpretation of the above considered partial solutions of equations (EP-(G$ 



First, we stress that at M = 0, i.e. in the Minkowski spacetime, equation (f23j) coincides 
with (j21j) and we obtain 

o mr sin 9 

m 3 = — . (46) 

<->2 

Naturally, in this case the restriction on the sin 9 does not appear, because the geometric 
point r = and the plane 9 = 71/ 2 are not connected with any physical source. 
Relationship (|4fi|) describes the known motions of the non-proper centers of mass of 
a spinning body. 

It is important that expression (J3H|) which was obtained for M^O does not pass into 
f!46|) at M — *■ 0. This situation is similar to the case for the equatorial circular orbits 
in a Schwarzschild field which was investigated in jKJj (relationship (}46|) at = 7r/2 
corresponds to equation (j2*Uj) from [10]). Therefore, by the arguments discussed in [10] 
we conclude that expression (JHHj) describes the motion of the proper center of mass of a 
spinning particle. 

In this connection we point out the fact of importance. At r = 15M/7, when by 
(JIHJ) 9 = 7r/2, according to equation (|2*2*|) we have S\ = 0, i.e. spin is orthogonal to 
the equatorial plane of the particle motion. The expression for the particle's orbital 
velocity for the equatorial circular orbits with 2M < r < 3M in a Schwarzschild field 
was written in [10! , equation (24). It is easy to check that this expression coincides with 
the corresponding expression following from equation at r = 15M/7, namely 

, 5mM 

Using (JHJ), (JHIJ), (|1D)) it is easy to see that for the considered non-equatorial 

circular orbits the relationships (w 4 ) 2 ^> 1, (ru 3 ) 2 ^> 1 take place. That is, the velocity 
of a spinning particle on these orbits is ultrarelativistic, as well as on the equatorial 
circular orbits which were studied in [40J. 



3.5. The case of the shortened Mathisson-Papapetrou equations 

In some papers, devoted to the investigations of the Mathisson-Papapetrou equations, 
instead of strict equations (0) their shortened form was used |H1EI|, namely 

m Ts uX = -\ U " SP(TR ^- (47) 
It means that condition (j2J) is imposed on the Mathisson-Papapetrou equations. As we 
pointed out in section 2, the replacement of equations (JTJ) by equations (|4Tjl eliminates 
the superfluous solutions of the strict Mathisson-Papapetrou equations. It is important 
in this approach that condition (J7j) must be sutisfied. 

In the context of our investigations it is interesting to verify (i) do equations (J2J), (}4Tj) 
under condition (j^J) have the non-equatorial circular solutions, and (ii) does condition 
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(J7J) is satisfied on such possible solutions. Therefore, in this subsection we shall consider 
equations (g7|), ©, © with relationships (HUJ), (HU), fl]©, (JUj), 1|22J . 
Two non-trivial equations of set (|4Tj) with A = 1 and A = 2 are 

3M 



3M 



-m 3 ^ sin 6 1 , 



From f49|) we have 

3 3M5i 



u S\ sin 9. 



u 



mr 3 cos 6 1 



By (HHJ), flUD, (El 



m z r° 



5M 



r 



1 I tan 2 9. 

r 



Substituting (joHj) . (joT|) in condition (j3T?j) we write 

18M/S'i\ 2 A 3M 
mr 



1 



tan 2 6> = 1. 



r \mr J \ r 

According to (|52j) it is necessary r > 3M. 

Using ©-(HU), (EDI), (ED we find the relationship between Si and S Q 



q2 _ c2 
°0 — °1 



1 - 



3M 



sin 2 0+1 



2M 



cos 2 fl 



X 



x 1- 



2M 



-2 



By JSl), (ESD 



o2 



m 2 r 2 9M 



(m 4 )- 2 cos- 2 ^ 



, 3MV . 2 /> A 2M \ 
1 ) sm 2 6>+ l ) cos 2 9 



x 



(48) 
(49) 

(50) 

(51) 

(52) 



(53) 



x 1 



2M 



5M\ 
2 sin" 2 9. 



r J \ , ■ (54) 

Relationship (|34*j) shows that condition © is satisfied if and only if the value r is close 
to 3M and cos 2 # is close to 0. Then equation (|54jl can be written as 



C2 

m?r 2 



361 + % 



(55) 



where 



0\ = 1 , 02 = COSb', 



< 5x < 1, |5 2 | < 1. 



(56) 
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By (J56j) condition (J52)) takes the form 

o2 r2 

4i = ( 57 ) 

Taking into account (|H7j) we obtain from (jHOjl . (f5Tj) 

So, shortened Mathisson-Papapetrou equations (|4"7|). © at condition Q have the 
non-equatorial circular solutions with the non-zero components of the particle velocity 
determined by |J5%]) . According to (JSHJ), (JSHJ) the relationships (« 4 ) 2 3> 1, (ru 3 ) 2 ;§> 1 
hold, i.e. the velocity of a spin particle on these orbit is ultrarelativistic. 

By the direct calculations it is not difficult to check that the partial circular solutions 
of shortened Mathisson-Papapetrou equations (J2J), ©, from this subsection satisfy 
condition (jJJ) due to (jHEj). 

3.6. The value of the expression S^ V P V on the non- equatorial circular orbits in a 
Schwarzschild field 

All partial solutions of the Mathisson-Papapetrou equations considered in this section 
were obtained under supplementary condition (j3J). At the same time the explicit 
expressions above allow to estimate the value of the left-hand side of condition (J3J 
on these solutions. 

Using equations (J22j) . (J22|l . (|34"j). it is easy to obtain the relationships 

S 3 ^ = 0, = 0. (59) 

Other two expressions from the left-hand side of condition (@J), namely S lu P u and 
S 2u P u , are not equal to 0. However, these expressions are proportional to e (where e is 
determined by (jHJ). So, for the small e the value {S^P^l is much less than 1. This fact 
may be usefull for further investigations of the solutions of the Mathisson-Papapetrou 
equations in a Schwarzschild field under condition (J3J. 



3. 7. Conclusions to section 3 

Thus, the strict Mathisson-Papapetrou equations under the Mathisson-Pirani 
supplementary condition have the partial solutions describing the non-equatorial circular 
orbits of a spinning test particle with constant latitude in a Schwarzschild field. The 
region of existence of these orbits by the radial coordinate is determined by expression 
([41)1. or approximately by (}4"2*|) . The dependence of the angle 6 on r is determined by 
expression (}36|) . or approximately by (jIHJ). According to (|43|) . for a spinning test particle 
hanging under the equatorial plane in a Schwarzschild field the maximum value of the 
corresponding angle 7r/2 — 9 is equal to ~ 47° and this value is achieved at r ^ 2.35M. 

The shortened Mathisson-Papapetrou equations have in a Schwarzschild field the 
solutions describing the non-equatorial circular orbits as well. However, the region of 
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existence of such orbits is much less than in the case of the strict equations. Namely, 
by ()56|) r is in the small neighborhood of r = 3M, and it/ 2 — 9 is close to 0. 

All above considered non-equatorial orbits of a test spinning particle in a 
Schwarzschild field are highly relativistic. 



4. The second-order Mathisson-Papapetrou equations for non-circular 
equatorial motions in the Schwarzschild coordinates r, (p 



The equatorial ultrarelativistic circular orbits of a spinning particle in a Schwarzschild 
field have been considered in [JU]. Because any small perturbation deviates the orbit 
of a spinning particle from the fixed circular orbit, it is interesting to investigate the 
possible non-circular equatorial motions of such a particle in a Schwarzschild field. Here 
we focus our attention on same aspects of this question. 

According to equation (10) from [10], for any equatorial motions when spin is 
orthogonal to the motion plane the relationship 



5*2 = ru^So 



(60) 



holds (S\ = 0, S3 = 0). Taking into account (J60j) . we write three nontrivial equations 
(HJ as 



S r 1 



2M 



(uW - it s u 4 ) + u 1 



m — 3Squ 3 u 4 I 1 



/ 2M\ x 4 . 3 3S M x 3 . 4 3M / M 

— 3S 1 ) u u u H u u u 1 

\ r J r r \ r 



3M 



r V r 



x 1 



(u 1 ) 2 ^ 4 + rS ( 1 



2M 



3M 



/ 3\3 4 ^0-^ 
[U ) u — — 



2M 



1 - — ^ uV) 3 +mTL, «V+?^ ( 1 - — ) S «V = 0, (61) 



r 



So, 4-1 -4 ix 6S M / 2M 
— (u tt — u u ) — I 1 



r 

3S n M 



±v± .4 

u 



2M 



-1 



-1 



„, 3 ,s n (1-^1 M V 



3S 



2M 



2M 



4\2 



| 6MS ^ _ ^ 2.1/ 



r / 



l/„,3\2„,4 



li (ti ) U 



2S n M / 3M 



^ M 1 (M 4 ) 3 +2mr 3 3 M 1 M 3 = 0, (62) 



S„r ( 1 - ™V' («W - S V) -3S„ (l -H)(l- ™'l „V,V - - - 



r / 



3S m 3 [ 



(1--) («': 



. 4/ 3S M 3 4 , 3MS / MV, 
+ « 4 (m + — - — A 4 + — ^ (1 I II 



r / 



2M 



-3 
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-rS 1 



2M 



x U V(u 4 ) 2 + 



3M 



3S n M 



mV)3 _3SoM 



2M 



3M 



2M 



-i 



« V + 2mr*« 1 M 4 = 0. 



4 ,.1„.4 
14' 



(63) 



Using the condition = 1 it is easy to check that among three equations (|BT |l -(|B ^l 
there are only two independent equations. 

Further transformations of equations ()61 )) -()63 )) are connected with the consideration 
of the known first integrals of the Mathisson-Papapetrou equations in a Schwarzschild 
field, the energy E and the angular momentum L of a test particle [S] 

DS^ 1 

E = mui + guUfj, + -gu,iS , 



-mu 3 - g 33 u^ 



DS 3lM 1 
as 2 



By (|oTH) expressions (jo^|) can be written as 

2M 



E = m 1 



+ w 3 r>V - 2r? 3 (u 1 )V] 

2M 



L = mr u 



2„.3 



+2r 4 4 ( M 1 )V- M 4 r>V 3 ] 

From expressions (}65|) . (}66j) we have 



(64) 



(65) 



(66) 



M 3 U 3 



U ' 



^fT + 7T r >^- 2r ?3 



3 «V 



E 



M u' J in a 



u u 
4 / 2M 



r 2 m 1 rSV> m 1 



1 - 



r 



« 4 = «< + ^r^V 3 -r>V 3 



1 1 1 a/3 



rSou 1 



(67) 



mr m 3 1m 4 / 2M 
+ — — r 1 

Here we put u l ^ (the case of circular orbits was under investigation in [40J). 
After the differentiation of equations (JBTj) . with respect to s we obtain 



(68) 



u^u 1 - iiV 



1 



u 1 



u 



a/3 L 



E 
rSn 



M o m 

2 + -5" 

r z rbo 



1 - 



2M 



11 



Some partial solutions of Mathisson-Papapetrou equations in a Schwarzs child field 15 



x 



2u 3 (I* - Its) 



E 



rS (uiy 
E Mu 3 



2 [it V + w 3 r^tiV 
2ikf 



X 



L mru 3 1 
rSo So r 



2M 



it 



T rV 4 - ^r^j^uV 



Wr 3 AlM W 



1 



2M 



it 



V l 



4M 



14 2M 14 
w it —it it 



+ 



M 1 rS'o 
mr« 3 1 



2M 



-w 



u 



a/3 1 



Taking into account equations (|57 jl -(|55 jl . we write equation (J53j) in the form 



m 1 - -fM 1 ) 2 -2r I 1 



3M 



f« 3 ) 2 - - I 1 



3M 



(69) 



rE 

Of) 



2M 



u 4 = 0. 



(70) 



It is easy to check that equation by (fBTf). ((BSj) coincides with equation (fTOj). From 
equation (J5T|) according to (J57j) . (|55|) we obtain 



u 3 + 



i s 
u u 



/ 3M\ (w 3 ) 2 £ fi 2 , 3N21 



-- ( 1 

r 



3M\ w 



(M + Lu 3 ) 1 U 4 = 0. 

nil. qrot J- .S„ \ r- / 



r J u 1 ru S\ 



2M 



(71) 



Using the condition w^w^ = 1 we write equations (|7D|). (|7T|) in the coordinates r and tp: 



r = h 2r ( 1 

r 



3M 



2 r # 



+ -(1 

r 



3M 



+ 



rSo 



r 2 + 1 



+ r 1 



2M 



3M\ 3 £ 



r J f rSof 



1 + rV) 



r(H-rV) 



1/2 



(72) 



+ 



m + 
rSof 



f 2 + 1 



2M 



(1 + rV) 



1/2 



(73) 



Thus, all possible equatorial world lines and trajectories (except the circular orbits) of 
a spinning particle in a Schwarzschild field are described by equations (|72|). ()73|). 
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We note that according to ()64|) at So — ► we have approximately E = mu^, 
L = —mu 3 , i.e. the expressions for the geodesic motions. Substituting these values in 
dZ2J>9 f73j) it is easy to see that the coefficients at the terms in both these equations 
which are proportional to l/S become equal to 0. That is, at S — > the right-hand 
sides in f)72jl and (J73J) are finite. Besides, it is not difficult to check that if the linear 
by spin corrections in expressions ()64|) are taken into account, then the right-hand sides 
of equations (J?2J), (|7H)) coincide with the corresponding expressions following from the 
equatorial geodesic equations. 



5. The concrete cases of highly relativistic essentially non-geodesic 
non-circular motions 

Here we shall confine ourselves to the cases of the non-circular motions when the initial 
value of the particle's radial velocity is much less then the initial value of its tangential 
velocity. Such a choice is determined by the known fact that the gravitational spin-orbit 
interaction decreases for the growing r as r -3 . That is, just on the circular or close to 
circular orbits a spinning particle feels the maximal effect of this interaction. 

We also point out that the ultrarelativistic circular orbits from [30] are chosen as 
the basic orbits in this section, i.e. the dynamics of the deviation of a spinning particle 
from these orbits is under investigation. 



5.1. The equations in the non-dimensional quantities 



For further calculations, instead of equations (|T2)) . (fTS)) which are written in the variables 
r(s) and ip(s) we shall use their representation through the non-dimensional quantities 



T 



M ds ds 

Then from equations (J72J), (f73|) we have 



P 



M 



dY 
dr 



Y< 



+ P 1 



v 

+ - 

P 



Y' 2 + ( 1 - - 

P 

YZ 



dZ 
dr 

1 



1 + Z 2 p 2 ] 
Z 2 + 1/p 2 



2Z 2 

1/2 



pZp 



p 



+ p 



Y 



^ = Y, 
dr 

^ = Z 
dr 



Z p 

p 

1/2 



(74) 



(75) 



(76) 

(77) 
(78) 
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where 

ME L , , 

t*=-Q-, V= ^- ( 79 ) 

Do Oq 

The small quantity e in (J76j) is determined by (jSjl. In the following we shall put e > 0, 
without any loss in generality Then according to (|50|) S 2 > 0. 

In the terms of quantities (|74)1. the condition that the initial radial 4- velocity is 
much less by the absolute value than the tangential 4-velocity can be written as 

\Y \<^p \Z \. (80) 

Different values of the parameters \x and v in equations (J75)) . (|76j) at the fixed 
initial values Po,fo, Y and Z correspond to the solutions describing motions of different 
centers of mass. It is important to pick out /i, v corresponding just to the proper center 
of mass. This question we shall consider below. 

5.2. The choice of the values fi and v for ro = 3M 

First, we put that the initial value of the tangential velocity is equal to the velocity of 
the spinning particle on the circular orbit of radius r = 3M considered in jlU]. Then 



(81) 

Using equations (p>4"j) and notations (J7§J), by the direct calculations we can obtain the 
values /i c and v c for the circular motion with r = 3M : 

+ ^ + (82) 

^-3W/*(l-f e + i^). (83) 

According to the analysis from [IHj, relations ()81|) - (|%3*|) correspond to the motion of the 
proper center of mass of a spinning particle. 

Naturally, in the cases of the non-circular motions, when Yq 7^ 0, the values /1, 
v describing the motion of the proper center of mass differ from f!82j) . (|H3*|) . However, 
because we are restricted by conditions (J80|) . l|8T|) . i.e. because we shall consider only 
the orbits close to the fixed circular orbit, we put 

/i = kxpc, v = k 2 v c (84) 

where the coefficients k\ and k 2 are close to 1. For the approximate estimates of values 
ki, k 2 we propose such a way. Let us consider the expressions for dY/dr and dZ/dr at 
the initial moment r = 0. From equations (J73J), (|75jl by ([50 )1 -(|S1 |) we obtain 



dr K ' 2 e 2 



/, , h ! 1 + 



55) 
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■ < 86 > 

According to notations (ffij) . expressions (J85|) . (}86|) give the initial values of the radial 
and tangential acceleration components. It is easy to see that the right-hand sides 
of equations ([85)1 . (JHUj) essentially depend on kx,k 2 and e. It is important that the 
right-hand side of equation (J85|) contains the term with the large quantity e~ 2 . In this 
connection we remember the corresponding expressions following from the Mathisson- 
Papapetrou equations in the Minkowski spacetime. Namely, it is known that in this 
case just the expressions for the accelerations of the non-proper centers of mass are 
proportional to (So/mr)~ 2 , in contrast to the acceleration of the proper center of mass 
which is equal to (see, e.g., |1BI EDI; the value So/mr in the Minkowski spacetime 
is analogous to e from (|SJ) in a Schwarzschild spacetime). Here we can suppose that 
the similar property posses the solutions of the Mathisson-Papapetrou equations in a 
Schwarzschild spacetime. Namely, that the term in (jH^J) which is proportional to e~ 2 
determines the motions of the non-proper centers of mass, whereas for the description 
of the proper center of mass it is necessary to put 

fci - (l + h « 0, (87) 

i.e. to decrease the effect of the term with e~ 2 on the initial value of the radial 
acceleration. 

We stress that the reflection above concerning the choice of //, u for the proper 
center of mass is not strict (its role is heuristic) and must be verified in the concrete 
computer calculations. We shall do it in the next subsection. 

It is easy to see that under condition (jS7|) the effect of the large quantity e~ 5 / 2 on 
the tangential acceleration in equation (|HB|) is decreased as well. Besides at the condition 

h « 1 (88) 

the influence of the large quantity e~ 3 ^ 2 in (fKHj) is decreased. Because on the circular 
orbit dY/dr = 0, dZ/dr = 0, it is naturally to suppose that at condition (J5UJ) expressions 
(jHSJ), (fHOj) cannot be too large in the case of the proper center of mass. 

Therefore, according to relationships (|H7|) . (JHSj) . in computer integrations of 
equations (ffH|) - (ffH|) we start with the expressions 

h = 1, k 2 = I 1 + ^^e) . (89) 



5.3. The results of the computer integration of equations J7<5| )-f7£|) at r = 3M 

Typical results of the integration of equations (|75|) - (|75|) at (jSlj) . (}%9"j) are summarized 
in figures 1, 2 in terms of the radial velocity dr/ds = Y, the value r/M = p and the 
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angle if at s = 10~ 4 . By (|81|) for this value e we have Z m —44. In figures 1, 2 the 
interval of integration of equations ()75 |l -()78 |l by r is [0-0.2]. According to figure 1 in 
this interval the angle ip is changed from to —8 radians that corresponds to ~ 1.3 
revolution around the Schwarzschild mass. The plot of v?(t) is n °t presented in figure 
2 because it practically coincides with the corresponding plot in figure 1. The plots 
of Z{t) is absent in figures 1, 2 because in the considered interval of r the value Z is 
practically constant. 

Figure 3 shows the solutions of equations (J7T)j) - (J7H|) in the case when the first relation 
from (jSHJ) is violated, namely at k\ = 1 + 2 X 10~ 7 , whereas the second relation from (|89|) 
remains the same. The oscillation motion of the non-proper center of mass is visible 
in figure 3. More generally, the computer integration of equations (ff3|) - (fTSj) shows that 
the similar situation takes place for other concrete values k±, hi which violate condition 
(jHHJ). In this connection we conclude that relation ()89)) are correct for the description of 
the motion just of the proper center of mass. 

In subsection 5.7 we shall discuss the results of subsection 5.3. 
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Figure 1. Y = dr/ds,p = r/M and (f vs. r by the strict Mathisson-Papapetrou 
(MP) equations for e = 1CT 4 , p = 3, Y Q = 2.5 x 1CT 3 . 
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Figure 2. Y and p vs. r by the strict MP equations for e = 10 4 , po = 3, Yq = 10 




0.002 0.004 x 0.006 0.008 0.01 




0.002 0.004 x 0.006 0.008 0.01 



Figure 3. The examples of the oscillation solutions of Y and p by the strict MP 
equations for e = 10~ 4 , p = 3, Y Q = 2.5 x 10~ 3 . 



5-4- The results of the computer integration of equations J75l )- J7^ ) at r = 2.5M 

It is pointed out in [3U] that in a Schwarzschild field the ultrarelativistic circular orbits 
of the proper center of mass of the spinning particle are admissible at 2M < r < 3M 
as well. In contrast to the circular orbit of radius r = 3M, which follows both from the 
strict and shortened Mathisson-Papapetrou equations, the orbits with r < 3M follow 
from the strict Mathisson-Papapetrou equations only. Similarly to the above considered 
deviation of a spinning particle from the orbit r = 3M, it is possible to study the 
deviation of a particle from the orbit with r < 3M. The initial equations are the same, 
namely (j75j) - (f7Sj) . 

As a typical case, let us consider ultrarelativistic motions which deviate from the 
circular orbit of radius r = 2.5M. By the results of [IHj, we found for this circular orbit 

Z = _ 2 V2 5 -i/4 £ -i/2 1 1 _ W5 £ _ a0 23 £ 2 ] ; ( 90) 
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(91) 



//, = 7 • 5- 3 /V2V 1/2 ( 1 - 



70 



Vr 



55/ Wi + S 



V2 



250 



(92) 



(the corresponding expressions for the circular orbit with p = 3 were written in (|8ip- 
©). 

By the procedure similar to the described above (see equations (|Sl jl -(|55 jl ). we can 
obtain the expressions for fi, v corresponding to the proper center of mass when Yq ^ 
at condition (JHUj) : 

100, 
821 



H = /j, c 



I - — r 2 (l + 87.5Y 2 )- 1 



(93) 



z/ = v r 



l_^2y 2(l + 87.5^-1 



821 



(94) 



Figure 4 shows typical results of the integration of equations fl73 j) - (|7H|) at conditions 
(|93|) - (j94|) for p = 2.5. These results are discussed in subsection 5.7. 




Figure 4. Y, p and ip vs. r by the strict MP equations 
100. 



for e = 10- 10 ,p = 2.5, Y a = 



Some partial solutions of Mathisson-Papapetrou equations in a Schwarzs child field 22 
5.5. The case of the shortened Mathisson-Papapetrou equations 

For the equatorial motions of a spinning particle in a Schwarzschild field equations (j4~ 
can be written in the coordinates r, ip as 



3M\ M -WIS 



o . 

— 2 V 



r 2 +(l 



2M 



r = ip r I 1 



ip = — rip. 
r 

At So = these equations coincide with the geodesic equations. 
In quantities (fTIj) equations (jHSJ), (j^)) can be written as 



nl/2 



(1 + rV) 



(95) 
(96) 



- 3e— 

P 2 

dz 
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dr 
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YZ 




P 





lV2 



2 „2^ 



1 + z> 



dip 
dr 



Z. 



(97) 
(98) 
(99) 
(100) 



An example of the integration of equations (|9*7 j) - ()100|) is presented in figure 5 . 
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Figure 5. Y and p and vs. r by the shortened MP equations for e = 10 4 ,po 
3,F = 2.5 x lCT 3 . 



5.5. T/ie case of the geodesic equations 

In notations (JUJ) the geodesic equations coincide with equations (|9*Tj) - (jl00|) when e 
in dHZj) is equal to 0. It is interesting to compare the corresponding solutions of 
the Mathisson-Papapetrou and geodesic equations. Figures 6-8 we shall use for the 
discussion of the results following from figures 1, 2, 4. 
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5.7. Discussion 

First, let us compare figures 1, 5 and 6 which show the corresponding solutions of 
the strict Mathisson-Papapetrou equations lj75j) - (j7B]) . the shortened equations (|H7j) - 
(jlUU|) and the geodesic equations respectively, with the same initial values of po,Y , Z : 
Po — 3, Yq — 2.5 x 10 -3 , Zq = —44, at e = 10~ 4 . As figures 1 and 5 show, on the r-scale 
of order 0.2 the corresponding plots of Y and p are closer. (The plot of <p(t) is n °t 
presented in figure 5 because it is not distinguished from <p(r) in figure 1). Whereas the 
plots Y(t) and p(r) in figure 6 essentially differ from the corresponding plots in figure 
1: according to figure 1 the spinning particle moves away from a Schwarzschild mass, 
while by figure 6 the spinless particle falls onto a Schwarzschild mass. We stress that the 
difference of the radial coordinate of both the particles Ar becomes comparable with r 
after less than 2 revolutions around a Schwarzschild mass. 

The initial values in figures 2 and 7 are: po = 3, Yq = 1CT 2 , Z = —44, at e = 1CT 4 . 
These figures show the situation when both the spinning particle and the spinless particle 
move away from the Schwarzschild mass but with different velocities. 

Let us compare figures 4 and 8. Here e = 10~ 10 ,p — 2.5, Y = 100 and, by 
expression Zq = —9.5 x 10 4 . According to figure 4 for the r-scale of order 

5 x 10~ 5 the spinning particle moves away from a Schwarzschild mass, reaches the 
value approximately p = 3 and after that moves toward a Schwarzschild mass. Whereas 
according to figure 8 the spinless particle monotonously falls onto a Schwarzschild mass 
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and for the r = 1.2 x 10~ 5 reaches the horizon surface p = 2. We point out that the 
plots of <p(r) i n figures 4 and 8 show that these processes correspond to less than one 
particle's revolution around a Schwarzschild mass. 

6. Conclusions 

The analytical solutions of the Mathisson-Papapetrou equations from section 3 and the 
numerical solutions from section 5 have a common feature: they describe such a type 
of ultrarelativistic motions of a spinning particle when the interaction of spin with the 
gravitational field acts as a repulsive force. It is connected with the fact that in all above 
considered cases the signs of Sg and dip/ds are opposite (the same relation holds for the 
circular orbit with r = 3M [40J). Naturally, when the signs of Sg and dip/ds coincide the 
spin-gravity interaction acts as a attractive force (this type of ultrarelativistic motions 
we shall consider in another publication). 

It is important that both the strict Mathisson-Papapetrou equations at the 
Mathisson-Pirani supplementary condition and the shortened variant of these equations 
admit the partial solutions in the form of the ultrarelativistic non-equatorial circular 
orbits with constant latitude in a Schwarzschild field (however, the region of existence 
of these orbits in the second case is much smaller, see equations (|5l)])). Moreover, 
on these solutions the expression \S^ U P U \ in the Tulczyjew-Dixon condition is much less 
than 1 due to the small value e which is determined by (JBJ). 

Naturally, the ultrarelativistic circular orbits of a spinning particle in a 
Schwarzschild field are unstable. The dynamics of the deviation of a spinning particle 
from the circular ultrarelativistic orbit with r = 3M in a Schwarzschild field caused by 
the non-zero initial value of the radial velocity is investigated in section 5. 

The result of importance is that spin can influence the shape of a trajectory 
very strongly, as compare to geodesic trajectory, for the short time of the particle's 
motion. Namely, for less than one or two revolutions of a spinning particle around a 
Schwarzschild mass the difference of the radial coordinate Ar of the spinning and spinless 
particle becomes comparable with the initial radial coordinate of these particles (see, 
for example, figures 1 and 6, or 4 and 8). We stress that this effect of the considerable 
space separation of particles with different spin takes place just for the small value e at 
the ultrarelativistic velocity with the Lorentz 7-factor of order e -1 / 2 . 

In jTJj we read: " The simple act of endowing a black hole with angular momentum 
has led to an unexpected richness of possible physical phenomena. It seems appropriate 
to ask whether endowing the test body with intrinsic spin might not also lead to 
surprises". Now we can answer this question in the positive sense, at least in the 
theoretical plane. 

Perhaps, the further analysis of the gravitational ultrarelativistic spin-orbit 
interaction, as well as of the gravitational ultrarelativistic spin-spin interaction (for 
example, in a Kerr spacetime), will be useful for more fine investigations of the 
gravitational collapse and cosmology problems. 



Some partial solutions of Mathisson-Papapetrou equations in a Schwarzs child field 26 



References 

Mathisson M 1937 Acta Phys. Pol. 6 163 
Papapetrou A 1951 Proc. R. Soc. A 209 248 

Horvathy P A 2003 Acta Phys. Pol. B 34 2611 (Preprint hep-th/0303099) 



Corinaldesi E and Papapetrou A 1951 Proc. R. Soc. A 209 259 

Das A 1957 Progr. Theor. Phys. 17 373 

Micoulaut R 1967 Z. Phys. 206 394 

Lawrence J K 1969 Acta Phys. Austriaca 30 313 

Mashhoon B 1971 J. Math. Phys. 12 1075 

Wald R 1972 Phys. Rev. D 6 406 

Rasband S 1973 Phys. Rev. Lett. 30 111 

Mashhoon B 1975 Ann.Phys. 89 254 

Tod K, de Felice F and Calvani M 1976 Nuovo Cim. B 34 365 

Epikhin E N and Mitskievich N V 1976 Acta. Phys. Polon. B 7 543 

Bovyn M P 1976 J. Phys. A 9 245 

Hojman R and Hojman S 1977 Phys. Rev. D 15 2724 

Carmeli M, Charach Ch and Kaye M 1977 Phys. Rev. D 15 1501 

Abramowicz M and Calvani M 1979 Mon. Not. R. Astron. Soc. 189 621 

Fuchs H 1983 Ann. Phys. 40 47 

Plyatsko R M 1988 Manifestations of the gravitational ultrarelativistic spin-orbit interaction (Kyiv: 

Naukova Dumka; in Ukrainian) 
Svirskas K, Pyragas K and Lozdicne A 1988 Astrophys. Space Science 149 
Yee K and Bander M 1993 Phys. Rev. D 48 2797 
Mino Y, Shibata M and Tanaka T 1996 Phys. Rev. D 53 622 
Tanaka T, Mino Y, Sasaki M and Shibata M and 1996 Phys. Rev. D 54 3762 
Suzuki S and Maeda K 1997 Phys. Rev. D 55 4848 
Suzuki S and Maeda K 1998 Phys. Rev. D 58 023005 
Saijo M, Maeda K, Shibata M and Mino Y 1998 Phys. Rev. D 58 064005 
Plyatsko R 1998 Phys. Rev. D 58 084031 
Semerak O 1999 Mon. Not. R. Astron. Soc. 308 863 
Mohseni M and Sepangi H R 2000 Class. Quantum Grav. 17 4615 
Suzuki S and Maeda K 2000 Phys. Rev. D 61 024005 
Bini D, Gemeli G and Ruffini R 2000 Phys. Rev. D 61 064013 
Mohseni M, Tucker R W and Wang C 2001 Class. Quantum Grav. D 18 3007 
Plyatsko R and Bilaniuk O 2001 Class. Quantum Grav. 18 5187 
Stuchlik Z and Hledik S 2001 Phys. Rev. D 64 104016 
Hartl M D 2003 Phys. Rev. D 67 024005 

Hartl M D 2003 Phys. Rev. D 67 104023 
Barrabes C and Hogan P 2004 Class. Quantum Grav. D 21 405 
Bini D, de Felice F and Geralico A 2004 Class. Quantum Grav. D 21 5427 

Bini D, de Felice F and Geralico A 2004 Class. Quantum Grav. D 21 5441 

Bini D, Cherubini Ch, Geralico A and Mashhoon B 2004 Preprint gr-qc/0411098 
Karpov O B 2004 Preprint gr-qc/0406002] 
Kao J-K and Cho H T 2005 Phys. Lett. A 336 159 
Plyatsko R 2005 Class. Quantum Grav. 22 1545 



Chicone C, Mashhoon B and Punsly B 2005 Preprint gr-qc/0504146 

Pryce M 1948 Proc. Roy. Soc. A 195 62 

Moller C 1949 Commun. Dublin Inst. Advan. Studies A 5 3 

Moller C 1949 Ann. Inst. Henri Poincare 11 251 

Moller C 1972 The Theory of Relativity (Oxford: University Press) 



Some partial solutions of Mathisson-Papapetrou equations in a Schwarzs child field 



[44] Fleming G 1965 Phys. Rev. B 137 188 
[45] Pirani F 1956 Acta Phys. Pol. 15 389 
[46] Tulczyjew W 1959 Acta Phys. Pol. 18 393 
[47] Dixon W G 1964 Nuovo Cim. 34 317 

Dixon W G 1970 Proc. R. Soc. A 314 499 

Dixon W G 1970 Proc. R. Soc. A 319 509 

Dixon W G 1973 Gen. Rel. Grav. 4 193 

Dixon W G 1974 Phil. Trans. Roy. Soc. 277 59 
[48] Frenkel J 1926 Z. Phys. 37 243 
[49] Mathisson M 1937 Acta Phys. Pol. 6 218 
[50] Weyssenhoff J and Raabe A 1947 Acta. Phys. Pol. 9 7 

[51] Synge J 1960 Relativity: The General Theory (Amsterdam: North-Holland) 

[52] Taub H 1964 J. Math. Phys. 5 112 

[53] Ragusa S and Bailyn M 1995 Gen. Rel. Grav. 4 193 

[54] Bonnor W 1977 J. Phys. A 10 1673 

[55] De Felice F 1979 Phys. Lett. A 69 307 



